function f = mle_reflection_rsc(v,l,g)
s_l = v(1);
n_l = v(2);
o_l = v(3);
e_l = v(4);
d_l = v(5);
s_g = v(6);
n_g = v(7);
o_g = v(8);
e_g = v(9);
d_g = v(10);
i = [12 11 10 9 8 7 6 5 4 3 2 1 0];
% The log-likelihood function is rescaled to avoid convergence issues with
% step size.
f = 10^5*(-1*log(prod((factorial(12)./(factorial(i).*factorial(12-i)).*(s_l.*(1-e_l).^i.*e_l.^(12-i) + n_l.*0.5.^i.*0.5.^(12-i) + o_l.*d_l.^i.*(1 - d_l).^(12-i))).^l))-1*log(prod((factorial(12)./(factorial(i).*factorial(12-i)).*(s_g.*(1-e_g).^i.*e_g.^(12-i) + n_g.*0.5.^i.*0.5.^(12-i) + o_g.*d_g.^i.*(1 - d_g).^(12-i))).^g)));
end